JOURNAL OF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 67D, No. 3, May -June 1963 

On the Statistical Theory of Electromagnetic 
Waves in a Fluctuating Medium (I) 

K. Furutsu 

Contribution from the Central Radio Propagation Laboratory, National Bureau of Standards, 

Boulder, Colo. 

(Received December 7, 1962) 

The subject of electromagnetic wave scattering by a randomly varying medium is 
reviewed giving special emphasis to the technical method of approach. The symbolic 
representation of Maxwell's equations is introduced to make it easier to survey the whole 
subject and to formulate the equations. The Feynman diagram method is applied to the 
computation of the correlation of the fields at different points in space to any order of 
approximation. The differential equation to be satisfied by the latter correlation function 
is also derived from another point ot view. Then 1 he theory is developed on the "renormali- 
zation" of the constants, i.e., the effective propagation constant in a fluctuating medium 
and the effective coupling constant between the field and the medium, etc.; the explicit 
expression of the former is obtained to the first order of approximation. The. dispersion 
relation is derived as a connected problem. In Part II of this series of papers, a fundamental 
theory of statistics of the electromagnetic field in a, fluctuating medium will be developed 
In Part III, a few applications to tropospheric scattering will be given. 

1. Introduction 

The entire subject of electromagnetic propagation through a fluctuating medium has been 
investigated theoretically and experimentally by many authors. Even though the Born 
approximation might be valid for the scattered electromagnetic waves in many eases, there 
still remain several theoretical problems to be solved such as the effect of earth diffraction, 
the effect, if any, of using very narrow beam antennas, the 4 restrictions for the validity of 
simple formulas for the scattered power, etc. As to the theory of multiple scattering, there 
are also many fundamental problems to be solved. 

The purposes of this series of papers are first to remark thai an almost perfect corre- 
spondence exists between the statistical treatments of the electromagnetic field in a fluctuating 
medium, which could be anisotropic, and the treatments in quantum field theory, and thus 
that the many powerful methods used in the latter, e.g., Feynman's diagram method, the 
method of renormalization, the met hods used in solving bound-state problems, etc., are available 
when solving the statistical problems in the former (Part I). Then a fundamental statistical 
theory of the field in a fluctuating medium will be developed (Part II); the field and the medium 
are inevitably quantized there, and the reason for keeping a high correspondence between the 
statistical theory and the quantum field theory will naturally be understood. 

The contents of part I can be summarized as follows: In section 2, a symbolic representa- 
tion of Maxwell's equations is introduced [Furutsu, 1952, 1956; Marcuvitz, 1962]. This 
representation is very convenient for studying the subject from a broader point of view, and 
for achieving a unified and brief equation formulation. In section 3, the theory of the power 
transmitted to the receiver through an arbitrary medium, including anisotropic media, is 
extended. Here the antenna patterns of the transmitter and receiver are represented in 
terms of the respective equivalent electromagnetic current distributions on suitable surfaces 
enclosing the transmitter and the receiver. These distributions could refer directly to pure 
electric currents of all the antenna elements which may be distributed on a surface of some 
extension, or they could be derived from the electromagnetic fields themselves, as occurring, 
e.g., on the surface of the mouth of a parabolic antenna. In the same way, the cross product 
of the amplitudes of two different receivers at arbitrary points in space is also derived in terms 
of the antenna patterns of the receiver and transmitter. Jn section 4, the expression for any 
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cross correlation of fields is derived with the aid of the results of section 3. All the equations 
are first transformed into their Fourier representations, and next Feynman's diagram method 
[Feynman, 1949a, 1949b; Dyson, 1949a, 1949b] is extensively used to any order of approxima- 
tion with respect to the correlation function of the refractive index. This diagram method 
has been one of the most convenient techniques in quantum electrodynamics. Here the same 
technique is introduced and, in section 5, the discussions are completely devoted to the "renor- 
malization," i.e., the effective medium constants in a fluctuating medium, the effective coupling 
constant between the field and the fluctuating medium, etc. 

In section 6, another approach to the multiple scattering problem is shown. It is used to 
derive directly the differential equation for the correlation function between the fields at two 
arbitrary points in space. The equation takes a simple form, and there are many possibilities 
to solve it. 

2. Symbolic Representation of Maxwell's Equations 

In order to survey the subject from a broader point of view and also to facilitate the 
analytic representations in the following sections, it is convenient to represent the six com- 
ponents of the electromagnetic field in a unified form. In this section, Maxwell's equations 
are expressed in symbolic form, and the Green function is properly defined for later 
convenience. 

Here we employ the following notations: Latin subscripts assume values ranging from 1 
to 3, and a repeated index is to be so summed. The coordinate vector in space is denoted by 
x t = (x) . 

Using the conventional notations, Maxwell's equations are given in Gaussian units (time 
factor e io)t ) bv 

rot H-i(u/c)&E=4TrI/c, rot E+i(u/c)nH=0. (2.1) 



(2.2) 
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the operator "rot" as a three-rowed matrix is expressed by 

d 



rot= Sibi=(sd), dj=i 



dxj 



0=1,2,3). 



Thus, (2.1) can be expressed in the form 

/0 1 



(*d) 





which takes the form 
Here 



[Pi(*d)-J#=j. 



/4tt 



and i= 



(2.3) 

(2.4) 

(2.5) 
(2.6a) 



are one-column matrices of 6 elements, and 



£=g{(i+P3)s+(i-P 3 M 



(2.6b) 



a 6X6 matrix and, in the, case of isotropic media, the elements 8 and n of the latter, as well 
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as those of the two-rowed Pauli's mat rices 

"-(',' } -(I "o) -C -D' (2 - 7) 

represent 3X3 unit diagonal matrices multiplied by the indicated quantity. 

On the other hand, in the case of anisotropic media, 8 and jjl are generally 3X3 matrices, 
which can be represented as functions of s x , s 2j and s 3 . Therefore S and n in general do Jiot 
commute with the s/s (see (2.10)). 1 

Equation (2.5) finally takes the form 

h£i-m=3 } Ji^PiSi. (2.8) 

The matrices s t and p, are all Hermitian and especially 

yt=il=-7i, (2.9) 

where 7* and y] are the complex conjugate and transposed matrices of y u respectively. Equa- 
tion (2.9) is evident from the explicit matrix expressions of (2.2). 

The matrices s t and p L satisfy the following algebraic relations which are often very con- 
venient to use: 

s i s j s k +s k s j s i =d ij s k +8 jk s i , x i Xj—Xj« i = ir 1 e ijk s k , 

4 = s 2 i + 4+«l=2, pi=pi=p3=l, piPi>=-- —f>2pi=ipz, etc. (2.10) 

Here e ijk is the antisymmetrical tensor corresponding to ei23=l; its elements are -f 1 or — ] 
according as the number of permutations necessary to change the order of subscripts ?', j, k, 
into the order 1, 2, 3 is even or odd. From (2.10), it especially follows, for any two orthogonal 
vectors m and // in the ordinary three-dimensional space, that 

(nyf=n 2 (ny), (ny)(my) 2 = {m 2 —('my) 2 } (ny) , (ny)(my)(ny)=0, 

(ny) 2 (my) 2 = (my) 2 (ny) 2 = n 2 m 2 — (m X n> • y) 2 > (2-11) 

The traces of these matrices are frequently used in actual computations, and the following 
formulas hold (appendix 1): 

Tr'(s,)=0, Tv'(s i s j )=28 ih Tr'( W *)=* _1 eu*, Tr'(s i s j s k s l ) = 5 i Ai+M". (2.12) 

Here Tr' stands for the trace in the snbspace of s/s matrices. 

The Green theorem in this representation takes the following form: Let \p" and 1// be 
arbitrary continuous ^-vectors represented by a one-row and a one-column matrix of 6 elements, 
respectively, in the space 2, then with (dx)=dxidx 2 dxz 

f {r / l(y^)-kW~^ f [-(yb)-kW}(dx)= -i\ V'(ny)yt,'da. (2.13) 

Here the right side represents a surface integral over the surface a of S, and n=n(o) the inward 

space unit vector normal to a, and d operates on the coordinates (x) on its left side. The 
excitation of the electromagnetic field will be performed either directly by electric currents, 
or indirectly by the use of a device such as a waveguide antenna. In the latter case, it will, 
in practice, be excited by one predominant waveguide mode. However, these two ways of 
excitation involve no difference whatever in the representation of the field: 

1 For example, in an ionized gas of plasma frequency co„/2ir under the influence of static magnetic field H or corresponding gyrofrequency 
w* J r//2 7 r with collision frequency v [Furutsu, 1952], 

S-l-wJMw-tV+As)}]- 1 . 
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Let \{/(x,x') be the 6X6 tensorial Green function satisfying given conditions in space and 
at boundaries. We then get, combining (2.8) with a result derived in appendix 2, two different 
representations for the differential equation for this function, i.e., 

[(yd)-kMx, x') = *(x, x')[~ (yh -k]=5(x-x'). (2.14) 

Here the right side is to be interpreted as the coefficient in front of a 6X6 diagonal unit matrix. 
On the other hand, \p(x, x') constitutes a nondiagonal 6X6 matrix, the jth column of which 
satisfies Maxwell's equations in which the right side of the ^th equation is b(x— x') instead of 
zero. The symbol d' =i<)/dx / here refers to differentiation with respect to x\ 
Then, the solution ipix) of (2.8) is given by the one-column matrix 

Hx) = fHx,x')j(x')(dx'). (2.15) 

On the other hand, the field which is excited by a source inside the surface a can be expressed 
in terms of the tangential components of the field on a. Thus, applying (2.13) to the space 
outside a, while putting ^ = ^{x) (the solution to be obtained) and \l/" = \p(x', x) (the Green 
function), we find, taking into account the outward propagating condition for \p(x) at infinity, 
and also using (2.14), 

^{x f ) =i \p(x\ x) (ny)\f/(x)da ) 

J or 

or 

+(x) = f4>(x,x')Ux')(dz'). (2.16) 

Here 

Mx)=id(x,*)(m)+(x), (2-17) 

while 5(x, a) is defined as the one-dimensional <5-f unction of a variable which may be identified 
with the distance to a, that is a function different from zero only on the surface a. 2 It may 
be remarked that %«(x) depends only on the components of the electromagnetic field that are 
tangential to the surface a. 

Equations (2.15) and (2.16) take the same form, and thus £«,(#) may be interpreted as the 
effective electromagnetic current accounting for the field inside a. 

In this symbolic representation, the Poynting vector P t takes the form 

Pi=£-ry*P={- [Re (E*XH)]i- (2.18) 

87T 47T 

Here ^* stands for the one-row matrix representing the complex conjugate of \p. Decomposing 
k into its Hermit ian and anti-Hermitian parts, 

k = K-i(47r/c)a, (2.19) 

both k and a (not to be confused with the same symbol for the surface a) being Hermitian, 
we readily find by the use of (2.8) and (2.6) that 

div P + ^W=-Re (E*I). (2.20) 

Thus the anti-Hermitian part of k contributes to the waves influenced by the dissipative 
properties of the medium. 



? For an arbitrary function fix), it holds that 

ff(T)8(x, a) (dl)=f a f(x)da. 
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On the other hand, from (2.17), we derive, using the relation (ny)y i (ny)=n i (ny), 

f&sr&>(dx)=nW,* r ), (2.21) 

if ^ is normalized by 3 



1 



^*(//7)^/t7=l. (2.22) 



The function 5(cr", O is the delta function of the "distance" between two "parallel" 
surfaces <j" and a' (according to t he usual definitions of these quantities for neighboring 

surfaces). 

In view of (2.18), the normalization (2.22) implies that the total power radiated outwards 
from the surface a equals c/Sw power units. Here, the surfaces a and o" in (2.21) are assumed 
to be two arbitrary elements of a continuous set of surfaces, no element of which intersects 
another. 

For later convenience, we shall deduce the Green function S(x) for the case of a homo- 
geneous medium which could be anisotropic; then, from (2.14), 

[(yb)-k)S(x)=8(x), (2.23) 

with the condition of outward propagation at infinity. Let S(t) be the Fourier transform 4 
of S(x) as follows: 

/» 00 

S(t)= e Wx) S(x)(dx), (dx)=dx 1 dx 2 dx Z} (2.24) 

then, from (2.23), 

[(yt)-k]S(t) = l, S(t) = [(yt)-k]-\ lm(fc)<0. (2.25) 

Here, hn (i)<0 indicates that all the eigenvalues of k to be used in (2.23) should have an 
infinitesimal negative imaginary part. 

On the other hand, the adjoint Green's function S(x) is defined by 

S(x)=S(-z)=-±r, fe +t ^S(t)(dt), {<lt) = <lU<lUdt„ (2.26) 

which satisfies, with the condition of outward propagation at infinity, 

S(x)[-(yb)-k] = 6(x). (2.27) 

In the special case of a homogeneous isotropic medium, k commutes with all the s t matrices 
(but still does not with p,), and thus, using the algebraic relation (ydf = £> 2 (yd) (see (2.11)), 
which involves (yt) 3 —(yt)t 2 , we find 

[(yt)-k][(yt)+ {(^0 , -*H-i /, }i- 1 ]=* , -fc' , ', 

in which the new matrix k and the new scalar k' are defined by 

F 2 =M=(o>/c)V, k = Pl kp l = k\^- l>r (2.28) 

Hence, from (2.25), 

S(t)=g(t)(?-k'Y\ g(t) = (yt)+{(yty-?+k' 2 }k-\ (2.29) 



3 This normalization is possible only when the total radiated power away from the surface a into the direction n is positive. 
* The notation t will be used exclusively for the Fourier transformation variable, and, in this paper, any function /(0 will stand for the Fourier 
transform of the space coordinate function f(x). 
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The function S(f) has poles at t 2 =k' which can be avoided in further integrations by 
assuming also an infinitesimal negative imaginary part for k' '; the latter corresponds to the 
condition of outward propagation at infinity. Hence, from (2.24) and (2.29), we finally have 

S ^=^£j- u ' x) S(t)(dt)=g(b) J^^y-k , y i e- i ^\dt)=g{b)S^ «-«''}. (2.30) 

Here r=\x\ and bj=ib/dxj. 

3. Transmission Power Gains 

We consider here a transmitter system which is assumed to be in the space 2! enclosed 
by the surface <ri and, in the same way, a receiver system in the space 2 2 enclosed by the surface a 2 . 
The conductivity loss in space and in the systems is first neglected, involving a Hermitian k. 
Let yp\{x) and fcix) be the solutions satisfying the given boundary conditions in space and at 
the boundaries, whose fields on the surface a^ or a 2 are respectively given by ^(o-,-) and are 
normalized according to (2.22) 5 by (fig. 1) 



I 



#(nY)fc&r=l (t=l,2). 



(3.1) 



Thus, since T T ,=7*=— 7, and df= — b f , 

Kyd)-k]Mx)=*t(fi)[-(y*)-k]=o, (3.2) 

the latter heing the complex conjugate of the former. However, the transposed expression 
of the second representation yields, 

[(Yd)-fr]«(aO=0. 

Hence, ^* constitute no solutions of Maxwell's equations unless k=k T . This holds for an 
isotropic medium ; the x//* then are exactly the solutions satisfying the given boundary conditions 
except for the propagation directions which are just inverse to those of \p t . 

Although (3.2) is sufficient in the case of isotropic media, it is necessary to introduce the 
adjoint equation of (3.2) in order to include the case of anisotropic media, i.e., 



*<[-(7&)-*] = 0, [(7&)-*Mf = (*=1,2), 



with the normalization 



- f fi(ny)tfd<r= f r,(ny)tid<r=l. 



(3.3) 



(3.4) 



Here ^ have to satisfy the same boundary conditions and propagation condition at infinity as 
yp t . It is to be noticed that ^* are the set of solutions of Maxwell's equations whose propagation 
directions are just inverse to those of ^. In the special case of k = k T , we have ^<=^«, and the 
normalization (3.4) agrees with (3.1). 



5 In this paper, (<n) will stand for the coordinates of arbitrary point bounded on the surface <n. 




n 2 n, 




Figure 1. Illustrations of the surfaces <ri, o->, oV and a? 
for (3.1) and (3.6). 
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Now suppose that both the transmitter and receiver use waveguides in each of which 
only one propagational mode can exist, and that the antenna systems are so well matched 
that there are no reflected waves at all for these modes. Any cross sections of these waveguides 
will be denoted by a[ and a'>, respectively, corresponding to the surfaces <ri and a 2 enclosed by 
them; the predominant modes will be marked \f/i(a[) and ^2(02) inside the respective waveguides, 
\f/i{x) and \l/ 2 (x) being their continuations outside. 

The problem is, then, to determine the field strength or the power transmitted into the 
receiver waveguide by the excitation of the mode wave in the transmitter normalized by (3.1). 
Putting \l/ / = \p l (jt) and \p" = \l/ 2 (x) in (2.13), and identifying the space 2 with the whole space 
outside the surfaces o^ and a 2 , we have, on taking into account that both \p x and \f/ 2 are solutions 
of (3.2) and (3.3) that represent waves propagating outwards at infinity, respectively, 






f 2 (ny)ypid<r=0. (3.5) 



Here we could replace the surfaces 0-1 and a 2 by a[ and 0-3, respectively. 

On the other hand, deep in the waveguides, say at <r' 2 and a[, we know ^(cr,) and foWi) to be 
waves entering into the second and the first waveguides, respectively; therefore 

fcW^sitffM), &(eri)=MfW), etc. (3.6) 

Here a n and b ]2 are numerical constants. Multiplying by fa((r 2 )(n 2 y) on the left side of the 
first equation of (3.6), and by (7^7)^1 (<r[) on the right side of the second equation, and inte- 
grating both sides over a 2 and a[, respectively, we find, using also the normalization conditions 
of (3.1) and (3.4), 

js 2 (>iy)\p ] <fo=—a 2h ,^2(^7)^1^=612, (3.7) 

J o 2 J°i 

and, on account of (3.5), 



1 



f h(ny)xp { (Ia=— a 2l + b 12 =0, or a 2l = b V2 . (3.S) 

2+*l 

This agrees with the theory of reciprocity in the case of isotropic media. However, when re- 
ferring to (3.6), we should bear in mind that \j/ 2 (x) is not a solution of Maxwell's equations but 
of the adjoint equation (3.3). 

Using the expression (2.18) for the Poynting vector, and also the normalization condition 
(3.1) (involving a transmitted power of c/Stt units), the power W n transmitted into the surface 
02 per unit power of the transmitter takes the form 

W 2l =- f tfforJlM*, (3.9) 

which, on deforming the surface <r 2 into a 2 , yields after the substitution of (3.6) for \p { (a 2 ) 

W 2l = -\a 2l \ 2 ( ^ny)^p!a=\a 2l \ 2 , (3.10) 



or, in view of (3.7), 



W 2l = f fa(ny)\px<l<j 

\J<To 



(3.11) 



In the same way, in the inverse case of radiation away from the surface cr 2 , the power W x2 
transmitted into the surface &i takes the form 

W 12 =\ [ Mny)xp 2 dcr\ (3.12) 
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When the medium is isotropic, we have \f/ i =\f/ ij and thus, from (3.7) and (3.8), W 2 i = Wi 2 , 
as it should. 

From (2.16), 

Hence, we have finally, by deforming a 2 to a 2 , 

W 2l =\jjl 2 (x)*(x, x*)$ 9l (x'Kdx)(dx') \ 

Here, in accordance with (2.17) and (2.21), 

h(x) = —i8(x } a)^(x)(m) = i8(x,a)(ny)^(x) , 



(3.13) 



with 



^eor,yil>(dx)=n j b{<T' , ,a r 



(3.14) 



which is derived from (3.4). 

If we would like to have the cross product of the amplitudes of two different receivers 
labeled by the subscripts 2 and 3, the corresponding product W^-, ii=«2i^3i would be given, 
instead of (3.13), by 

'^^^'''L/T** 2 ^^^' »')fc 1 (»0(rfx)(rfeo]x[jy&,(x)^(a:, a5')& ri (x / )(d«)(da: / )]* (3.15) 

4. Fourier Transform and Feynman Diagram Method 

We assume here that the medium is changing randomly in space and time, but still has 
some correlation between two points in space which depends on space-independent averages. 
Also the time variation is assumed to be very slow as compared with the wave frequency. Let 
k=k J r 8k, where, for the time being, we put k =(k), the time average value of k and thus 
(3k) = 0. In an anisotropic medium, the k generally becomes a matrix operating on \j/. 

Then, from (2.14) fixing the Green function \k(x, x'), we find the equation 

[(yb)—k ]$(x f x')=8k(x)t(x, x')+8(x-x') k=k +8k, (4.1) 

whose solution is, using the Green function S(x— x f ) of (2.23) for the homogeneous medium k , 



iKx, x')= fs(x-x'')8k(x'0t(x'',x')(dx")+S(x-x')=Mx, aO+ftfo *)+fa(x f x') 



+ . . . (4.2) 



Here 



Mx, x')=S(x-x'), h(x, x')= f \dx")S{x-x")8k(x")S(x"-x'), 

fa(x, x f )= Udx") Udx ff/ )S(x-x f/ )8k(x / ')S(x /f -x f/ ')8k(x /,/ )S(x ,/ '-x / ), etc. (4.3) 

On the other hand, according to (3.13), the time average (W 2 i) of the transmission power 
gain W 2 \ is given by 

(W 21 )=jjjj(dx 2 )(d Xl )(dx^ (dxiXS^zWixl, x£g 2 (x£Xl 2 (x 2 )Mx2, x 1 )^ 1 (x 1 ))=J2T ij . 

(4.4) 

Here \p f is the Hermitian conjugate of \p an d defined as \k f (x[, x 2 )=comiplex conjugate of the 
transposed function yp T (x 2 , x[). Further, in view of (4.3), 
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r«= fjjj(dx 2 )(dx l )(d4)(dx' 1 )^* 1 (x' l ^Uxi, zl)& a (.zl)xl a (.Xt)M**, *i)M*i» (4.5) 
or, in an abbreviated representation, 

with 

Generally, 3k will have the form 

dk(x)=pAe(x). (4.6a) 

Here /3 is a 6X6 nondiagonal matrix in the case of an anisotropic medium, having the dimension 
of a propagation constant, and A£(x) is always a dimensionless scalar such as the fluctuating 
part of the dielectric constant of the medium. In the case of an isotropic medium, /3 has the 
form, in view of (2.6), 

r 3=i(l+ P3 )(co/c). (4.61)) 

Substituting the expression (4.6a) in (4.5), we see that the average value appears only 

through the term 

D(x h x 2 , x z , . . ., x n ) = (Ae(x 1 )Ae(x 2 ) . . . A£(x fl )), (4.7) 

which may be represented with the aid of its Fourier transform I)(ti, t 2j . . ., t n ), in the form 

00 

JJ. . . f(dt ;1 )(dt 2 ) . . . (dOe-' {,rt+, « + - +f » x » , Z?(*i, t 2t . . ., t n ). (4.8) 

On the other hand, the function I){.r u x 2) . . ., x n ) has two important properties: its 
invariance with regard to parallel translations and simultaneous rotations of the points x u 
x 2 , . . ., ./'„, and also the symmetry with regard to these points. Owing to these invariances, 
the functional form of U(t\, t 2 , . . ., t n ) is very restricted: from the translation invariance 
for Xi—^Xt+l, I being an arbitrary vector, we see by the use of (4.8) that the factor exp [— i(t\ 
+/ 2 + • • • -\-tn)l] should effectively be equal to 1 when multiplied by D(t u t 2 , . . ., t n ). 
This implies that, also taking into account the rotational invariance, the latter takes the form 

D(t 1} U, • • ., *n)=5(*i+* 2 + • . . +tn)D n (Q,U i ), (4.9) 

and thus the vector sum U+t 2 + . . . +t n should always be equal to zero. Here D n (t 2 iy t£ 3 ) 
stands for a function which is symmetric with regard to the indices i and j of the indicated 
independent variables. 

Substituting the expression (4.9) into (4.8), we finally arrive at 

D(x u x 2 , . . ., x n )= JJ. . . Udh)(dt 2 ) . . . (dt n )8(t 1 +t 2 + . . . +t H ) 

L °° XD n (t 2 u t i t j )e- Ht i x i +t 2 x 2+--+W. (4.10) 

Thus we also infer the invariance for coordinate inversion 

U ( X[ } x 2j . . ., x n )=LJ(xi, x 2 , . . ., x n ) . 
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In the same way, we define the Fourier transforms %,(() of £ a (x), an d £„(?) of £<,(#) by 

£„(*)= (2tt)- 3/2 f " f.(t)« + "' Jd (<ft); (4.11a) 

»/ — oo 

hence, if k T =k (see sec. 3) 

F.(0=^-*),* , =*. (4.Hb) 

The normalization condition (2.21) then takes the form 

P 8><t)yMtXdt)=nX<r" , »')■ (4-12) 

«y — oo 

Also, from (2.24), the Green function S(x) and its Hermitian conjugate function S f (x) are found 
to be 

S(z) = (27r)- 3 r S(t)e-*™(dt), S\x) = (2t)-* T SWe-^idt), S\t)=-S(t)\ k ^, (4.13) 

J — co J —co 

where lc* is Hermitian conjugate to k. 

It would be instructive to know the Fourier transformed form of T n . From (4.3) and 
(4.5), we derive, with the aid of (4.6a) and (4.7), 

xh 2 (x 2 )S(x 2 -x)pS(x-x 1 )^ 1 (x 1 ).D(x'--x). (4.14a) 
Here, according to (4.10), 

CO 

D(x'-x) = (AZ(x')A&(x))= N \dn(dt)8(t+t')D(t 2 )e- Ut '*'+ m == f " D{t 2 )e' mx, ~ x) {dt). (4.14b) 
Evaluating (4.14a) with the aid of (4.13) and (4.11a), we arrive at 

rn=JJJJJ(^D^ 

Xl^SiQffSit^^dih-to-^Dit 2 ). (4.15) 

Here, the terms 5(ti — t 2 — t) and 5 ( — t[-\-t' 2 -\-t) result from the integrations over x and x', 
respectively. 

Equation (4.15) finally takes the form 



in which 



T n =j^D(t 2 \f*(t)f(t)(dt) = (D) f , (4.16a) 

/» oo 

/(o= h^ti-^sih-^mti^ihWh) (4.i7) 



constitutes an ordinary scalar function. This integration will be evaluated in Part III for 
arbitrary functions for ^ and & . 

In the same way, for the cross correlation of the amplitudes at two receivers, we derive, 
starting from (3.15), 

W 32i u= ["_ D(t*)Mt)fn(t)(dt), (4.16b) 

which is the cross correlation in the first approximation with respect to D(t 2 ). Here j 2 \ is 
the function given in (4.17) and/ 3 i the corresponding function referring to the receiver labeled 
by the subscript 3 instead of 2. In Part III, we shall see that f*i(t)f 2 \{f) does oscillate more 
rapidly with / according as the distance \x z — x 2 \ increases, and thus the integral becomes smaller. 
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The integral (4.15) can be interpreted graphically in the following way (fig. 2a): On a 
sheet of paper, we first mark four points corresponding to £ ffl , | ff2 , £*- x and |j 2 , respectively, and 
the two points s and s' corresponding to the two scattering points. Then draw the solid line 
from the point £ ffl to s and from s to £<r 2 . Here the direction is taken away from the point £„,; 
it corresponds to the propagation with the wave-number vector t x away from £ ffl (connected 
with the transmitter, for instance) to the scattering point s. On the other hand, from the 
end point of this vector we draw a line towards the point £„ 2 ; it corresponds to the propagation 
into £ ff2 (connected with the receiver) with the wave-number vector t 2 . 

In the same way, we draw the solid line from the point |* 2 to $', and from s' to £* r But 
the directions of these complex conjugated lines are to be taken just inverse to those of the 
corresponding original lines. Hence, the new solid line is directed towards the point £* with 
the corresponding vector t[. FinahV we connect the points s and s' by a broken line with any 
prescribed direction of the vector t. 

Now there exists a one-to-one correspondence between the lines thus drawn in the diagram, 
and every factor in the integrand of (4.15). First of all, the factors ^(U) and £<r 2 (f 2 ) in the 
integrand fix the antenna patterns of the transmitter and receiver in the spaces enclosed by 
the surfaces <r x and <r 2 , respectively; they correspond to the points £ ffl and \ 2 in the diagram. 
In the same way, the factors £ x (t\) and & 2 (t' 2 ) correspond to the points £* x and |* 2 . Then the 
vector line U from the point £,, to s* in the diagram corresponds to the factor S(U) and, in the 
same way, the vector lines t 2 , t[ and t 2 to the factors S(t 2 ), S*(t[) and S f (t 2 ), respectively. The 
vertices s and s' are situated in accordance with the vector relations t=ti — t 2 = t[ — t 2 , thus 
expressing geometrically the effect of the momentum functions P8(ti — t 2 — t) and /3<5( — t[-\-t 2 J rt) } 
respectively. Finally the broken line represents the last factor D(f) of the integrand of (4.15). 

Since the mentioned factors are all 6X6 matrices or one-row or one-column vectors of 
six elements, their orders of succession in the integrand are important: from the right to the left 
these orders correspond toJ.be directions of the solid lines, thus beginning with the factor % a] 
and ending with the factor £ ff2 for the solid line considered first, while beginning with the factor 
|* 2 and ending with the factor gj^ for the complex conjugated second line. On the other hand. 
broken lines will always represent an ordinary scalar function, that is D(t 2 ) in the case under 
consideration. 

Though this is just an example of the way for representing the integrand of the multifold 
integral for T n by a corresponding diagram, its method is quite general and can be applied to 
all other cases. There is one such diagram for each T tj occurring in (4.4), and any of them 
comprises two solid lines, the one being the original line starting from £ ffl and ending at £ o2 , 
and the other the complex conjugate line starting from £j 2 and ending at £* ; in the case of 
the diagram for T ih there are j vertices on the original solid line from each of which one broken 
line is branching off; the corresponding property holds for the i vertices on the conjugate solid 
line. All these broken lines are connected to the "domain" of D which, including all of them, 
corresponds to the factor D(t h t 2 , . . ., tf+j) introduced in (4.9) (fig. 3). The 4 other corre- 
spondences of each solid line and each vertex in the diagram with the corresponding factors 
of the integrand are just the same^as explained above for T n . 
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Figure 2. Feynman diagrams for T u , To () , and T iy2 . 
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Figure 3. Feynman diagram for T31 in ter?ns of 

D(ti,t 8 ,t,,t4). 
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The number of broken lines or bonds of D is equal to the number of factors 8k involved in 
the integrand, and generally there are accordingly (? +j) bonds for the special term T fj . There- 
fore, for the second-order approximation terms, we have to do with the three terms T u , T 20 , 
and T 02 , the diagrams of which are given by figures 2a, 2b, and 2c, respectively. In the higher 
approximation of nth order, only the sum of terms T fj with i J rj=n would have a clear physical 
meaning, but in most cases, not the other individual terms. 

It remains to discuss the function D(x h x 2 , . . ., x n ) or DJfi. t t tj) which is of the nth 
order with respect to A£(x). Of course (A£(x))=0 and further, assuming a multivariate normal 
distribution of A6(x) values at different points, it is well known 6 that 

(Ae(^)A6(j 2 )Aefe) ■ . . Ae(x 2n )) = 2 p (AE(x 1 )Ae(x 2 ))(Ae(x z )Ae(x,)) . . . (A x 2n -. 1 )Ae(x 2n )), (4.18) 

(A£fe)Ae(x 2 ) . . . Ae(x 2n+1 ))=0. 

Here X p stands for the sum of all the terms obtained by the various combinations of different 
pairs composing the set of variables x u x 2 , . . ., x 2n . 

The right side of (4.18) is given in terms of the correlation function D{x) between only 
two points and this fact makes the computation very simple; first of all any diagram having 
an odd number of vertices yields no contribution to the result, and the function D(t u t 2 , . . ., 
U n ) can be divided into (2n— 1) (2n— 3) ... 1 independent parts depending only on the function 
D(f). 

For instance, the diagram of D(t u t 2 , t 2 , U) can be expressed as the sum of the three in- 
dependent diagrams of figures 4a, 4b, and 4c, and thus the diagram for T 3] is divided likewise 
into the three diagrams shown in figures 5a, 5b, and 5c. The contribution from each of these 
diagrams can be easily formulated in the integral form according to the correspondence rule 
stated above. 

The number of diagrams, however, rapidly increases with the order n, and the diagram 
method would become impractical. But, just as in the case of quantum electrodynamics, a 
renormalization method can be applied. The sum of contributions of all the self-energy parts, 
such as shown in figure 6a, and the sum of contributions of all the vertex parts, such as those 
of figure 6b, can then formally be taken into account by the use of an effective value of k, and 
another effective value of the matrix (3 introduced in (4.6a), respectively. This procedure is in 
accordance with the method developed by Dyson [1949a, 1949b]; we then only need to consider 
irreducible graphs, that is, graphs having no self-energy parts and no vertex parts at all. In 
the next section, the theory of renormalization, especially of k, will extensively be discussed. 
In particular, the imaginary part of the effective value of k which results from the renormaliza- 
tion, and which represents physically the attenuation of waves due to the scattering by the 
medium, will be computed to a first approximation, together with the corresponding dispersion 
relation. 



e See. for example, Middleton: Introduction to Statistical Communication Theory, ch. 7- sec. 7.3 (7.28), (McGraw-Hill Book Co.. 1960). 
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Figure 4. Decomposition of D(t,, t 2 , U, t 4 ) into the 
three independent diagrams in the case of a multi- 
variate Gaussian distribution. 



Figure 5. Three independent diagrams for T3 
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Figure 6 (a) Example of self-energy part; 
(b) example of vertex part. 
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5. Effective Medium Constant and Effective Fluctuating Medium-Field Coupling 

Constant (Theory of Renormalization) 

Here, (4.1) for the Green function is replaced by 

[(yty—keMx, x')=6k(x)f(x t x') + b(x-x'), k=k e +dk. (5.1) 

Hence, by averaging over the time, we obtain 

[(yV)-k § ](1,(2, x')) = (Sk(x)1,(x, x')) + 8(x-x'). (5.2) 

Here, k e is the time-independent effective value of k, which is different from k of (4.1), and 
which is to be defined by 

K(Hx 9 x')) = (yty(rKx, x'))-b{x-x'). (5.3a) 

Thus k e could be an operator operating on (\f/), but it will be assumed as a space- 
independent constant for a while. Since \^(./ // , x')) should depend only on the coordinate 
difference x n — x\ its Fourier integral representation will have the form 

(+(x'',x'))=^y 3 j W J e (t)e-"*''-*-\dt), (5.4) 

which is similar to the formula (4.13) for S(x). Here, according to (5.3a), S e (t) satisfies the 
equation 

l(yt) -k.]S.(t) = 1, S e (t) = [(yt) -keT l , (o.5) 

from which we conclude inversely that (\p(x, x')) satisfies the adjoint equation 

<*(*, x'))[-(yb')-k e ]=S(x-x'). (5.3b) 

It will be seen later that the k e in (5.5) is generally not constant but depends on t. This 
situation corresponds to replacing the left side of (5.3a) by 

k e (+(x, x'))->fk.(x-x")(iKx", x')){dx"), (5.6) 

where 

k e {x)= { ~^^k e {t)e-^{dt). (5.7) 

The Fourier transform of (5.3a) then leads to (5.5). 
Now, from (5.2) and (5.3a), we infer 

(5k(x)$(x, x'))=0, ($(x, x')8k(x'))=0, (6k(x)) = 6K9*0. (5.8) 

Here we put, as in section 4, 

(k)=k e + 5K = k . (5.9) 

By comparing (5.3a) with (2.23) we find, identifying the constant k in the latter equation with 
the present k 0J 

(Hx, x'))\ k ^ k =S(x-x'). 

315 



Equation (5.3a) is then solved for arbitrary k e by 

(xP(x, z')> = - fS(z-x")3K(W, x')){dx")+S{x-x'). (5.10) 

Here duty) should be interpreted as in (5.6). 

On the other hand, in view of the wave equations (5.3a) and (5.3b) for (\p(x, x')), the 
solution of (5.1) can be given in the form 

*(*, x') = Smx, z"))Bk(x>'Mz", x')(dx")+(Kx, *')> 

= fHx, x")5k(x")(t(x", x'))(dx") +(*(*, *')). (5.H) 

Here 6k(j) — 8k should be equal to 8lc(x) defined in (4.6a). Hence, 

8k(x)-8ic=0Ae (a?), (A£ (x)>=0, (5.12) 

and thus, substituting (5.12) in (5.8), 

teGKs, x'))+P(Ae (x)f(x, x'))=0. (5.13) 

It will be convenient to introduce now the "vertex function" Bfa; x z ) x 2 ), to be defined by 
(+(x h x 2 )Ae(x z )) = f(Hx h z'))B(z'i z'"\ z")(W, x 2 ))D{x'"-x z )(dx>){dx")(dx>"). (5.14) 
The Fourier transform of B will have the form 

B(x i; x 3 ; X2 )=^—f 8 (t 1 +t 2 +t d )B(-t 1} ««* {CW+(lrf+(W, (*i)(dW(*«)- (5.15) 

(2vr) b 

Here the 5-function appears because the left side of (5.15) should be invariant for a parallel 
displacement of the vectors x t , x 2j and x d , just as in (4.10). Applying the first expression of 
(5.11) to \l/(xx, x 2 ) on the left side of (5.14) and also using (5.12), we have 

(+(x u x 2 )At (x 3 ))=f(Hxi, x'))(8ktx')Hx', x 2 )Ae(x*))(dx') 
which gives, applying the second expression of (5.11) to yp{x' ', x 2 ) on the right side, 

fMzt, x'))$(Ae {x')Ae {&))&&, x 2 )){dx') 

+ f(+(x, x'))(5k(x')*(x', x")8k(x")Ae(x z ))(Kx", x 2 )){dx'){dx"). 

Thus we infer, since the above equation should be identical with the right side of (5.14), 

fB(x l ; x\ x 2 )D(x-x z )(dx)=p8(x 1 -x 2 )D(x 1 -x z )+(8k(x 1 )f(x u x 2 )8k(x 2 )AE (a*)). (5.16) 

Hence B becomes in the lowest order of approximation with respect to A£ 

B(xi; z 3 ; x 2 )^ph{x x — x 3 ) b (x 2 —x 3 ), 

which agrees with (5.15) if B(t u t 2 ) is given by the corresponding approximation 

£&,«<*,£ (5.17) 

When the values of A£(x) follow the Gaussian distribution, it will be proven in Part II 
of this series of papers that, for any functional /(A£[2]) of A£(x) over the space S 7 , 

(Ae(x)f(Ae[x]))= fr>(x-x') ^^/(a£[s])\ (&'). (5.18) 



7 The notation /(A<?[2]) means that the function / involves the variables A£(x) at all points of the space 2; for example, any waves are func- 
tional of the medium through which they propagate. 
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Here the functional derivative is defined as follows: When, for arbitrary variation <5A£(j), the 
change 8f can be represented in its first order by the form 



8 y(Ae [2 ])=J s {^/(A6[ S ])}aAe( 



x)(dx), (5.19) 

the factor of 8A6(x) in the integrand is defined to be the functional derivative with respect to 
AE(x). 

Applying the lemma (5.18) to the second term on the right side of (5.16), we find, taking 
into account (5.12) and also the relations 8 

8(8k(x))/8A£(x')=(38(x-x'), (5.20) 

H{x u x 2 )=j+(x u x')$bAt(;x')t{x', X 2 )(dx'), 
together with (5.8), that 

(8k(x^(x ly x 2 )8k(x 2 )Ae(x z )) = ^D(x z -x)(8k(x 1 )^(x u x)M{x, x 2 )8k(x 2 ))(dx). (5.21) 

Hence, it follows from (5.16) that 

B( Xl ; x,; x 2 ) =08(x 1 -x< i )8(x 2 -x 3 ) + (8k(x 1 )f(x 1 , x z )ffl(xt, x 2 )8k(x 2 )). (5.22) 

From (5.22), we can get B in any order of approximation by successive application of 
(5.11) and (4.18). 

Now, substituting the integral representations (5.4), (5.15), and (4.14b), and evaluating 
the integrations leading to impulse function, the integral (5.14) reduces to (fig. 7) 

(4,{x uX2 )AZ(x z ))=j^ C^« I ) J B(< 1 ,< 2 )^« 2 )Z?(« 1 -^)e- < («i (I i- I 3»+«2^3-2»)(d( 1 )(^2). (5.23) 
Hence, (5.13) becomes 

8k(*(x, x'))+ [[/3 fs § (ti)B(t lt * 2 )£((*i-f 2 ) 2 )(^i)] ~)3 S e (t 2 )e- i ^-^(dt 2 )=0. 

Thus, applying (5.6) and (5.7), we arrive at the Fourier transform of 8k, i.e., 

-dK(t)=fpS e (t 1 )B(t 1 , VDtth-ty)^). (5.24a) 

Here the Feynman diagram of 8k is represented in figure 8. 



8 The first equation of (5.20) follows from the identity 

AS(x) = \ 6(x-x')A£(.x')(dx'). 



■L 



|t,-t 2 



B(t„t 2 ) % S e (t,) >b 



\ t 2 \ D((t,-t 2 ) 2 ) 

Figure 7. Feynman diagram for {5.23). Figurf 8. Diagram for {5 .2 4a) . 
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In the same way, we derive from (5.10) the following Fourier transform S e (t) of 
(t(x",x')) (see (5.4)):"' 

S e (t)=S(t)-S({)8,i(t)S e (t)=S(t)[l-dK(t)S e (t)]. (5.24b) 

In a first order of approximation, S e (t) is given by 

S e ®~S(t)[l-dK(f)S(fi)], 

and, substituting (5.17) in (5.24a), <k(f) accordingly by 

In the case of an isotropic medium, the 8k({) defined in (5.24a) behaves as if it were a 
constant whenever it operates on the solution £(/) of the wave equation (compare (2.25)) 

[(yt)-kMt)=0, (5.26) 

which involves t 2 = k e k e = k' e . The reason for this behavior as a constant is the following: 
Since 8n(t) should be invariant with respect to space-coordinate rotations, it should be a 
function of t 2 , (yt), (yt) 2 , ((yt) 3 = t 2 (yt)), and thus it can be expressed in the form 

«Jc(0=M« + **i(f)[W)-tJ + ^(«»)[(7«)-ifcJ 2 . (5-27) 

Hence, when operating on %(t) of (5.26), 8k(() is equivalent to the constant of 8K Q (k' e ) which 
should be equal to k — k e according to (5.9). Thus 

8K (k' e 2 )=k -k e (5.28) 

gives the equation to determine k e . Also in the case of anisotropic media, we would have the 
same conclusion that <k(£) behaves as if it is a constant matrix whenever operating on the 
wave function of (5.26). 

5.1. Evaluation of the Effective Medium Constant for an Isotropic Medium 

(in First Order of Approximation) 

It will be interesting to compute the 8k as defined by (5.9), especially its imaginary part, 
in a first order of approximation in the case of an isotropic medium. Then, from (5.25), 
applying (2.29), 

-8K(t)=j™ (tl-tY^git^Dati-tyXdh). (5.29) 

Here, when the medium fluctuation is due to the dielectric constant, we can derive the following 
relations with the aid of (4.6b) and (2.29) : 

0=i(i+p 8 )(«AO, P 2 =(o>/c)p, ij8=(»Aoeft j8ff(Oi»=e- 1 {(*i7) a -«+fc / *}ft (5.30) 

and thus (5.29) takes the form 

faCO-ise- 1 p D((t 1 -m(dti)=-pe- 1 r (tf-i'V(7O 2 #((*i--0 f )(<«i), im(*')<o. 

%J — oo J — CO 

(5.31) 

Since the correlation function D(x) is an even real function of the space coordinates x, 
the Fourier-transformed function D(t 2 ) is a real function of t. Also fig(fi)P is real in so far as 
it is Hermitian, and thus has only real eigenvalues. 
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On the other hand, /•' and k' should have an infinitesimal negative imaginary part (see 
the end of see. 2). Hence, in the integrations, 

(tl-k ,2 )- 1 =P(t 2 1 -k ,2 )- 1 -iri8(tl-k /2 ). (5.32a) 

Here, P denotes the Cauchy principal part. Thus 

In. IMOH/^- 1 P d(t~-k' 2 )(yt ] YD((t l -t) 2 )(dt l ), 

%) — 00 

ReMoi-ise- 1 f " B{{h-midh)=-m-' \ " Pin-k^yKyhYDah-tmdh). (5.32b) 

*J — CO \) — 00 

The integral for Im [<5k(/)1 should be a function of t 2 , (yt) and (yt) 2 in view of the same in- 
variance discussed when dealing with (5.27); thus, when operating on the wave function £(/.) 
of (5.26), it behaves as a constant. In this case, the factor (yt x ) 2 of the integrand can be replaced 
by 

[r(0(7fe) 2 «0]/[f*«)W)] = 4«(l+cos»ff). (5.33) 

This result is easily obtained by using the method to he described in Part III. The right side 
of (5.33) constitutes the mean value for the two independent polarizations of £(£),and 6 is the 
angle between t t and t. Thus, putting (dti) = \ti\ 2 d\ti\dQ and dtt=2ir sin Odd, (5.32b) yields, 
when evaluating the integration with respect to |^|, 

Im (a*)=j8 \ k'*£- 1 f U+cos 2 d)D((2k' sin d/2) 2 )dQ. (5.34) 

Hence, from (5.28), we find 

Im (k e ) = -(3k/~ l £ {\ (1+ cos 2 d)a(d,k')dQ. (5.35) 

Here 

a(0, k , ) = (7r/2)k' 4 £- 2 B((2k' sin d/2) 2 ), (5.36) 

which is interpreted usually as the scattering cross section per unit volume from the view point 
of scattering loss. 

5.2. Dispersion Relation 

A well defined dispersion relation existing between 8k u the anti-Hermitian part of <k, and 
<k R , the Hermitian part, follows from (5.32b). This relation makes it possible to express <k R 
in terms of 8kj; by the change of variable ti—t=t', the right side of (5.31) takes the following 
form, when Im (k ; ) tends to zero: 

Equation (5.31)— pE" 1 [ " A^±D^j}l _ D(t' 2 )(dt'). (5.37) 

Here, e is an infinitesimal positive real number. From now on, the integral (5.37) will be 
understood as operating on £(f) of (5.26). Hence, the term { (t-\-f) • y} 2 in the integrand can be 
replaced by its mean value for the two independent polarizations of £(£), and thus by the right 
side of (5.33) with f 2 =k /2 . 

Since €>0, we recognize the integral (5.37) as an analytic function of k f in the lower half 
complex plane. Further, taking into account that 5k x , which is given by (5.34), will become 
the leading term of the left-hand member of (5.31) for k' -> + °°, we see that, for k / ^oo j this 
member tends to 

ijB(T/2)*' 8 e- 1 (D((k'd) 2 )di2 =ipT 2 k'e~ l f °° D(t 2 )tdt~k'. 

J k'->co JO 
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On the other hand, for k'—^0, (5.37) becomes 

z>(r 2 )(r/n+o[£']. 

J — oo 

Hence, in view of (5.31), it follows that 

(l*-$&~ l f~ D(P)(dt))/Jc' 2 ^8f(k' 2 ) (5.38) 

is a function of k' 2 which is analytic in the lower half plane and also on the positive real axis, 
while the singularity at the origin, if any, is integrable. Also, this function tends to zero if 
|t r | 2 -^oo in the mentioned domain. Hence 

^'^JTlS^ 0>-e>W), 
or, using (5.32a) with Im (&' )-»—€, 

5f(lc' 2 )=~-. (~ € P (s-k'Y^mds. (5:69) 

ITl J - m -ie 

Hence, taking the Hermitian parts on both sides of (5.39) in the limit e-^0, 

5/k(£' 2 H-- P Pis-k'Y'dUsyis, e->o. (5.40) 

7T Jo 

Here 8f R and 5/ T are the Hermitian and anti-Hermitian parts of 5/, respectively, and 5/i(s) is 
zero for s<0 because of (5.32b). By the substitution of (5.38) into (5.40), we obtain 

dK R (k ,2 )-l /36- 1 P D(t*)(dt)=~— P P(s-F 2 )- 1 ^ rfe. (5.41) 

o J -oo ttJo ^ 

Equation (5.41) is the dispersion relation which gives 8k r in terms of the known <k T - 
spectrum. 

Although this dispersion relation is derived for 8k in the first order of approximation, the 
corresponding formula would hold also for the exact 8k independent of any approximation. 

6. Differential Equation for the Cross-Correlation Function 

Another approach to the multiple scattering problem is as follows: Let us introduce the 
function G(x 2 , x 2 \xi, x[) defined by 

G(x 2 , x' 2 \x u x[)=(yl/(x2 J xi) • $(?&, xi)). (6.1) 

Here the sign * indicates the complex conjugate, and the product a • b the direct product of two 
matrices a and b. Hence, the function G is a 6X6 = 36 matrix, the elements of which equal 
the products of two corresponding elements of its matrix factors. Then, from (4.4), (W 2 i) 
is given by 

<w t 21 hJJJJi, 2 ^ (6.2) 

provided the vectors are suitably multiplied to G. 

Now, applying the wave equations (5.1) together with (2.3) and (2.9), we find 

[(y& 2 )-k e }'[(y<>2)-k*]G(x 2 , x 2 \x u x[)=(8k(x 2 ^(x 2 , Zi)-$k*(x0^*C4 x[)) 

+8(x 2 f -x' 1 )(8k(x 2 )rt/(x2 } x 1 )) + 8(x 2 -x l )(8k*(x^*(x2, x[)) + 8{x 2 -x x )8{x 2 -x[) . (6.3) 
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Here, according to (5.8) and (5.12), 



(6.4) 



If the medium fluctuations satisfy a multivariate normal distribution, we can use the lemma 
(5.18) and thus get for the first term on the right side of (6.3) 



+ \[Bk-P*D(i£—x)+P-8iPD(x2—z)] 



5&£(x) W^a^O'^W^O }/(&?) 



Here 



-/8-/S* fz>(x 2 — J?)Z?(Xa 



-a/)/ 



8AE(z)8AE(x') 

^^-^^(AeWAe^)). 



7t{*(^,Xi)-**W, »0} ){dx)(dx'). (6.5) 



As we have seen in section 5, /?(#) and 5k are of the second or higher order with respect 
AS. Hence, all the terms except the first on the right side of (6.5) are of the fourth order in 
AS, and thus may reasonably be neglected in most cases as compared with the first term which 
is of the second order. Hence 



(8k(x 2 )+(x 2) ^).6k*(/,)^(4, x[))^<P-P*D(x2-4)G(x 2) x' 2 \x u x[) } 



(6.6) 



and thus (6.3) becomes, remembering the vanishing of the second and third terms of its right 
side: 

[(7d 2 )-W-[(7d0-Wfe, aJk, x[)~P-p*D(x 2 -x' 2 )G(x 2 , £\x U 7$ +8(x 2 -x 1 )- 8(x' 2 -x[). (6.7) 

This constitutes the differential equation to be solved to get the function G. In the case of 
scalar waves, the corresponding equation reads: 



[A 2 + k 2 e ][A2+kf]G(x 2 , x 2 \x h x[) =cD(x 2 — x 2 )0(x 2 , x 2 \x u x' 1 )+8(x 2 —x 1 )-8(x2—x[). 



(6.8) 



The approximation used to derive (6.7) is equivalent to taking into account all the con- 
tributions from the graphs of ladder type like figure 9a but not those from graphs like figure 9b 

The investigation of the solution of (6.7) or (6.8) is also a powerful approach to the multiple 
scattering problem. 9 



9 The differential equation of this kind may he solved relatively easily hy the use of the relative coordinates u=X2— X\ and those of averages 

U=(X2+Xi)/2. 



X, Xj 

(a) 



>r 



X, X, 

(b) 



Figure 9(a) Diagram for ladder type interaction; 
(b) diagram for nonladder type interaction. 



7. Summary and Discussion 

In (his paper, the treatment of a statistical theory of electromagnetic fields in a fluctuating 
medium is primarily based on the graphical method introduced by Feynman [1949] for quantum 
electrodynamics. The latter also has been one of the most convenient methods in general 
quantum field theory. The theory is developed in parallel with that of Feynman as much as 
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possible, and a perfect correspondence between them is found in almost all kinds of problems 
occurring in our scattering theory, especially when the medium fluctuations follow the Gaussian 
distribution. For instance, the problem of getting the correlation between the fields at two 
different points in space (sec. 6) is equivalent to solving the Schrodinger equation for two 
mutually interacting particles, whose interaction term can be obtained in the same way as in 
quantum field theory [Salpeter and Bethe, 1951 ; Gell-Mann and Low, 1951]. The latter has a 
perfect correspondence with our theory especially when the medium fluctuations satisfy the 
above mentioned conditions. Many methods and approximations used in the quantum field 
theory [Levy, 1952a, 1952b; Klein, 1953] will therefore be available for solving the partial 
differential equation (6.7) for the most relevant correlation function; also the theory of 
renorinalization treated in section 5 has an almost perfect correspondence with that of the 
quantum field theory. 

The correspondence in question refers to the behavior of the electromagnetic field in a 
randomly changing medium, and to that of the quantized electron field which is in interaction 
with the quantized electromagnetic field; A£(x) behaves just like a quantized Boson field, 
especially when A£(x) has a Gaussian distribution so that (4.18) will hold. This might not be 
accidental, because the quantum theory states the problem and the answer only in terms of 
probability or statistics; any quantized field cannot be free from its own fluctuation. In Part 
II of this series of papers, the fundamental statistical theory for the electromagnetic field in a 
fluctuating media will be extended. The quantum-field theoretical treatment will inevitably 
be introduced there, and the basic relationship between the statistical theory of electromagnetic 
fields and the quantum field theory will be derived; many powerful methods used in the latter 
are expected to be available in the former in view of this relationship. In Part III, a few appli- 
cations to tropospheric scattering will be discussed, taking into account the boundary conditions. 



The author expresses his cordial thanks to Dr. H. Bremmer for his great interest and en- 
couragement extended to this paper. 

8. Appendix 1 

Introducing the matrices p i {i=\ ) 2, 3) by 

fr=2s?-l, /3 2 =2sf-l, ft=2«S-l, (Al) 

these prove to satisfy the relations 

(i^j) ftft-frft=0, 

(i=j), ftj8 2 3 =l, «=l,etc. 






Hence, applying the algebraic relations (2.10) and also the relation Tr(ab) = Tr(ba) for arbitrary 
matrices a and b, 

Tr'(s l )=Tr'(l3 1 s l )=Tr'(s 1 1 ) = hTr'(l3 l s 1 +s 1 l3 1 ) = O, 

Tr'(s 1 s 2 )=Tr'(fi 1 s l s 2 )=-Tr'(s 1 s 2 p 1 )=-Tr f (fi 1 s 1 s 2 )=-Tr'(s 1 s 2 )=0 l 

Tr'{s\)=\Tr'{s\+sl+si) = 2, 

Tr'(s i s j s k )=Tr'(—s k SjS i +8 i} s !c +& jk s i ) = —Tr'(s t s j s i )=2 Tr'[(s,Sj— SjS t )s k ]=-^ e m Tr'{siS k )=i- l e m . 

(A3) 

Here, the commutation relation in (2.10) is used. In the same way, since the trace is invariant 
for the matrix transposition and sf=—s ij 

Tr'(s i s j s k $ l )=Tr\s l $ k s j s i ) = ±Tr'[(s i s j s k +s k s j (A4) 
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Generally, usin^ the commutation relations for s t and fi t , the trace of any function of s L 
can be expressed in terms of the following traces: 

Tr'( Sl ) = Tr'(s ,(ij) = Tr'istsM = 7>'MA)=0 (i^jVfc), 

2V'(«^)=2« tf> 7>% W ,) =;-'.,;,,, 7V'(/3,) = + l, Tr'd3 t (3j)= -1 (tV./). 

9. Appendix 2 

The second equation of (2.14) is derived from the first one as follows. The Green's 
function \f/(x, x') is defined by the first equation 

[(yb)-m{x,x')=b{x-x') t 

which gives, operating [—(yd')-k] from the right side, 

[(t&) -*]*(», *0[-^ 

Hence 

[(yt»-k]{+(x,x')[-(yh-k}-8(x-x')}=0, 

which has the general solution 

H*, x')[-(yV)-k] = 5(x-x')+<t>(x, x'). 

Here <t>(x, x') is a solution of the homogeneous equation 

[(yd)-k]<t>(x } 0=0, 

which has to satisfy all given boundary conditions in space. Generally, there exists no such 
solution of a homogeneous equation. Thus (f>(Xj x')=0 and the second equation of (2.14) is 
obtained. 
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